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ABSTRACT 

We report the identification of cyelical changes in the orbital period of the echpsing 
dwarf nova Z Cha. We used times of mid-echpse collected from the literature and 
our new eclipse timings to construct an observed-minus-calculated diagram covering 
^ . 30 years of observations (1972-2002). The data present cyclical variations that can be 

' fitted by a linear plus sinusoidal function with period 28±2 yr and amplitude 1.0±0.2 

I minute. The statistical significance of this period by an F-test is larger than 99.9%. 

Tij" i The derived fractional period change, AP/P — 4.4 x 10""^, is comparable to that of 

' other short-period cataclysmic variables (CVs) , but is one order of magnitude smaller 

, than those of the long-period CVs. Separate fits to the first and second half of the data 

(N . lead to ephemerides with quite different cycle periods and amplitudes, indicating that 

the variation is not sinusoidal or, most probably, is not strictly periodic. The observed 
C~| . cyclical period change is possibly caused by a solar-type magnetic activity cycle in 

^ i' the secondary star. An incremental variation in the Roche lobe of the secondary star 

' I of ARl2/ Rl2 — 1.7 X 10""* is required in order to explain both the observed period 

^ • change and the modulation of the quiescent brightness previously reported by Ak, 

^ ■ Ozkan & Mattel. 

' Key words: accretion, accretion discs ~ stars: dwarf novae - stars: evolution - 

binaries: eclipsing - stars: individual: Z Cha. 



1 INTRODUCTION the continuous angular momentum loss are magnetic brak- 

„ „, . , ■ 1 , , „o , X ,. . ing via the secondary star's wind (for Port > 3 hr) and 

Z Cha IS a snort-period (Forb ~ 1-78 nr eclipsmg cata- , ,. ,r n ^ n t \ n^- irvoo\ 

, , .Z..... r , , , gravitational radiation (for Porb < 3 nr) (King 1988). At 

clysmic variable (CV). In these binaries, a late-type star i ^ • j i j_i i ^ i^ r n 

, , \ r-ni 1 n 1 r verv short periods, when the secondary star becomes luUy 

(the secondary) overfalls its Roche lobe and transfers matter , . / , ^ ^ ^ » ^ \ i i i . 

^ . , . , ^ , , . . T , degenerate (M2 kj 0.08 Mq), mass loss leads to an expan- 

to a companion wfaite dwart (the primary). In most CVs the ■ r • j_ 1 j_i 1 ^ 1 i^- 

^ \ f- J y gjQjj ^j-^jg g^g^j. g^jjjj j-everses the secular trend, resulting 

donor star has lower mass than the accreting star. Since con- , , r, ■ ■ ■ 1 1 • 1 tt ii 

° thereaiter m an increasing orbital period. However, the pre- 
servative mass transfer m such situations would lead to an 1 . . 1 , r . r, ■ • 1 • • • 1 

dieted mass transter rate alter this period minimum is low 
increase m the orbital separation (and therefore the cessa- 1^-12 -i\ 1 r r^xr ^ 1 ^ 

^ ^ (M2 ~ 10 Mq yr ) and tew CVs are expected to be 

tion of mass transfer via Roche lobe overflow), the existence , , . 1 1 i- ^ /nr i/ir,r\ 

, , , observed m such evolutionary stage (Warner 1995). 
of CVs as mass-transler binaries implies that they must con- 
tinuously loose angular momentum in order to sustain the 

mass transfer process. As a consequence, the binary should The secular evolution of the binary can in principle be 

evolve slowly towards shorter orbital periods (on time scales detected by measuring the changes in the orbital period of 

of 10* - 10^ yr). Possible mechanisms suggested for driving eclipsing CVs. Echpses provide a fiducial mark in time and 

can usually be used to determine the orbital period (and 
its derivative) with high precision. However, attempts to 

* email: bapQastro.ufsc.br (RB); chico@das.inpe.br (FJ); sim- measure the long-term orbital period decrease in CVs have 

biotica@starmedia.com (EO); sonja@mensa.ast.uct.ac.za (SV); been disappointing: none of the studied stars show a pos- 

pwoudt@artemisia.ast.uct.ac.za (PAW); msc@astro.keele.ac.uk itive detection of an orbital period decrease (e.g., Beuer- 

(MSC) mann & PakuU 1984). Instead, most of the well observed 
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Table 1. Log of the observations. 



Table 2. New eclipse timings. 



Date 


Filter 


Telescope 


Cycles 


Cycle 


HJD 


BJDD 


(0-C)| 


1995 Oct 23 


V 


1.6m/LNA 


130 874, 130 875 




(2,450,000+) 


(2,450,000+) 


(cycles) 


2000 Apr 14 


W 


0.6m/LNA 


152 817-152 820 


130 874 


0014.70294 


0014.70358 (4) 


+0.0115 


2002 Feb 12 


R 


1.9m/SAAO 


161795, 161797, 161799 


152 818 


1649.51443 


1649.51515 (4) 


-0.0043 


2002 Feb 13 


R 


1.9m/SAAO 


161810, 161812 


161 803 


2318.89047 


2318.89119 (4) 


-0.0072 



eclipsing CVs Qshow cyclical period changes (e.g., Bond & 
Freeth 1988; Warner 1988; Robinson, Shetrone & Africano 
1991; Baptista, Jablonski & Steiner 1992; Echevarria & Al- 
vares 1993; Wolf et al. 1993; Baptista et al. 1995; Baptista, 
Catalan & Costa 2000). The most promising explanation of 
this effect seems to be the existence of a solar-type (quasi- 
and/or multi-periodic) magnetic activity cycle in the sec- 
ondary star modulating the radius of its Roche lobe and, 
via gravitational coupling, the orbital period on time scales 
of the order of a decade (Applegate 1992; Richman, Apple- 
gate & Patterson 1994). The relatively large amplitude of 
these cyclical period changes probably contributes to mask 
the low amplitude, secular period decrease. 

Z Cha seemed to be a remarkable exception in this sce- 
nario. The eclipse timings analysis of Robinson et al. (1995) 
shows a conspicuous orbital period increase on a time scale 
of P/IPI — 2 X 10^ yr, not only at a much faster rate than 
predicted (~ 10^ yr) but also with the opposite sign to the 
expected period decrease. 

In this Letter we report new eclipse timings of Z Cha 
which indicate a clear reversal of the period increase ob- 
served by Robinson et al. (1995). The revised (O— C) dia- 
gram shows a cyclical period change similar to that observed 
in many other well studied eclipsing CVs. The observations 
and data analysis are presented in sectioiiQ and the results 
are discussed and summarized in section pT 



2 OBSERVATIONS AND DATA ANALYSIS 

Time-series of high-speed differential CCD photometry of 
Z Cha were obtained on October 1995 and April 2000 at the 
Laboratorio Nacional de Astroffsica, Brazil, and on February 
2002 at the South African Astronomical Observatory. The 
observations were performed under good sky conditions and 
covered a total of 11 eclipses while the target was in qui- 
escence. The time resolution ranged from 10 s to 30 s. The 
April 2000 run was made in white light (W). A summary 
of these observations is given in Table 111. Data reduction 
included bias subtraction, flat-field correction, cosmic rays 
removal, aperture photometry extraction and absolute flux 
cahbration. A more complete analysis of these data will be 
presented in a separate paper. Here we will concentrate on 
the measurement of the mid-eclipse times. 

Mid-eclipse times were measured from the mid-ingress 
and mid-egress times of the white dwaxf eclipse using the 
derivative technique described by Wood, Irwin & Pringle 
(1985). For a given observational season, all light curves 
were phase-folded according to a test ephemeris and sorted 
in phase to produce a combined hght curve with increased 



t i.e., those with well-sampled observed-minus-calculated (O— C) 
eclipse timings diagram covering more than a decade of 
observations. 



" Observed minus calculated times with respect to the linear 
ephemeris of Table ^ 



phase resolution. The combined light curve is smoothed with 
a median filter and its numerical derivative is calculated. A 
median-filtered version of the derivative curve is then an- 
alyzed by an algorithm which identifies the points of ex- 
trema (the mid-ingress / egress phases of the white dwarf) . 
The mid-eclipse phase, 4'o, is the mean of the two measured 
phases. Finally, we adopt a cycle number representative of 
the ensemble of light curves and compute the corresponding 
observed mid-eclipse time (HJD) for this cycle including the 
measured value of 0o. This yields a single, but robust mid- 
eclipse timing estimate from a sample of eclipse hght curves. 
These measurements have a typical accuracy of about 4 s. 

For Z Cha the difference between universal time (UT) 
and terrestrial dynamical time (TDT) scales amounts to 
19 s over the data set. The difference between the baricen- 
tric and the heliocentric corrections is smaller than 1 s as 
Z Cha is close to the ecliptic pole. The mid-eclipse timing 
have been calculated on the solar system baricentre dynam- 
ical time (BJDD), according to the code by Stumpff (1980). 
The terrestrial dynamical (TDT) and ephemeris (ET) time 
scales were assumed to form a contiguous scale for our pur- 
poses. The new eclipse timings are listed in Table ^ The 
corresponding uncertainties in the last digit are indicated in 
parenthesis. 

The data points were weighted by the inverse of the 
squares of the uncertainties in the mid-eclipse times. We 
arbitrarily adopted equal errors of 5 x 10~^ d for the op- 
tical timings in the literature (Cook & Warner 1984; Cook 
1985; Wood et al. 1986; van Amerongen, Kuulkers & van 
Paradijs 1990), and a smaller error of 2 X lO"'^ d for the 
timing of Robinson et al. (1995), half the error in our mea- 
surements, and combined the 74 timings to compute revised 
ephemerides for Z Cha. Table ^ presents the parameters 
of the best-fit linear, quadratic and linear plus sinusoidal 
ephemerides with their l-o formal errors quoted. We also 
list the root-mean-squares of the residuals and the value 
for each case, where V2 is the number of degrees of freedom. 

Fig. presents the (O— C) diagram with respect to the 
linear ephemeris in Table ^ van Teeseling (1997) reports 
that a positive offset of 90 s (plus an additional offset of 
0.0025 cycles) was needed in order to centre the x-ray eclipse 
of Z Cha around phase zero with respect to the ephemeris of 
Robinson et al. (1995). We used this information to assign a 
representative eclipse cycle to his measurement and to com- 
pute the (O— C) value with respect to the linear ephemeris 
in Table ^ His x-ray timing is plotted as a cross in Fig. ^ 
and is fully consistent with our results. 

The significance of adding additional terms to the lin- 
ear ephemeris was estimated by using the F-test, follow- 
ing the prescription of Pringle (1975). Not surprisingly, 
the quadratic ephemeris is no longer statistically signifi- 
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Table 3. Ephemeridcs of ZCha 

Linear ephemeris: 

BJDD = To + Po ■ -B 

To = 2440264.68070 (±4) d 



XZ 



■ 49.2, V2 



72 



Quadratic ephemeris: 

BJDD = To + Po ■ -B + c ■ 

To = 2440264.68088 (±8) d 
c = (+2.8 ± 1.1) X IQ-^* A 
xl, = 48.0, ,^2 = 71 



Po 
o"l ■ 

Po 

0-2 - 



-- 0.0744993048 (±5) d 
: 3.27 X 10-3 cycles 



= 0.074499300 (±2) d 
: 3.13 X 10"3 cycles 



Sinusoidal ephemeris: 

BJDD = To + Po • -E + A- cos [2-k{E - B)/C] 

To = 2440264.6817 (±1) d B = (120 ± 4) x 10^ cycles 
Po = 0.074499297 (±2) d C = (136 ± 7) x 10^ cycles 
A = (7.2 ± 1.0) X 10-* d (Ts = 1.09 x lO^^ cycles 



: 3.99, ^2 



10 



o 

X 



o 
o 



o 




20 40 60 80 100 120 140 160 
cycle number (x 1000) 

Figure 1. The (O— C) diagram of Z Cha with respect to the 
linear ephemeris of Table ^. The optical timings from the litera- 
ture are shown as open circles, the x-ray timing from van Teesel- 
ing (1997) is shown as a cross, and the new timings are indi- 
cated as solid circles. The dashed line in the upper panel depicts 
the quadratic ephemeris of Robinson et al (1995) while the solid 
line in the lower panel shows the best-fit linear plus sinusoidal 
ephemeris of Table ^. Best-fit linear plus sinusoidal ephemeridcs 
for the data on the first half of the time interval (B < 80 X 10'^ 
cycle, dotted curve) and for the last half of the time interval 
(E > 80 X 10^ cycle, dashed curve) are also shown in the lower 
panel. 



cant. From the ai and its values in Table |3| we obtain 
F(3, 69) = 130.7, which corresponds to a 99.99% confidence 
level for the linear plus sinusoidal ephemeris with respect to 
the linear fit. The best-fit linear plus sinusoidal ephemeris is 
shown as a solid line in the lower panel of Fig. y. 

However, the eclipse timings show systematic and sig- 
nificant deviations from the best-fit linear plus sinusoidal 
ephemeris. The fact that xi^ > 1 emphasizes that the linear 
plus sinusoidal ephemeris is not a complete description of 
the data, perhaps signalling that the period variation is not 
sinusoidal or not strictly periodic. We explored these possi- 
bilities by making separate fits to the first half of the data 
set (£ < 80 X 10^ cycle) and to the second half of the data 
set (is > 80 X 10'' cycle). The results are shown in the lower 
panel of Fig. |l| respectively, as dotted and dashed curves. 
The older timings point towards a longer cycle period with 
a lower amplitude (A= 1.13 min, Vcyc = 29 yr), whereas the 
more recent timings are best described by a shorter cycle pe- 
riod of larger amplitude (A= 1.23 min, Pcyc = 23 yr). These 
results indicate that the period changes are not sinusoidal 
or not periodic. 



3 DISCUSSION 

Our results reveal that the orbital period of Z Cha is no 
longer increasing as previously found by Robinson et al. 
(1995). Instead, the (O— C) diagram shows conspicuous 
cyclical, quasi-periodic changes of amplitude 1 min on a 
time-scale of about 28 yr. 

Cyclical orbital period changes are seen in many eclips- 
ing CVs (Warner 1995 and references therein). The cycle 
periods range from 4 yr in EX Dra (Baptista et al. 2000) to 
about 30 yr in UX UMa (Rubenstein, Patterson & Africano 
1991), whereas the amplitudes are in the range 0.1-2.5 min. 
Therefore, Z Cha fits nicely in the overall picture drawn from 
the observations of orbital period changes in CVs. 

If one is to seek for a common explanation for the cycli- 
cal period changes in CVs, then models involving apsidal 
motion or a third body in the system shall be discarded as 
these require that the orbital period change be strictly peri- 
odic, whereas the observations show that this is not the case 
(Richman et al. 1994 and references therein). We may also 
discard explanations involving angular momentum exchange 
in the binary, as cyclical exchange of rotational and orbital 
angular momentum (Smak 1972; Biermann & Hall 1973) re- 
quires discs with masses far greater than those deduced by 
direct observations, and the time-scales required to allow 
the spin-orbit coupling of a secondary of variable radius are 
much shorter than the tidal synchronization scales for these 
systems (~ 10* yr, see Applegate & Patterson 1987). 

The best current explanation for the observed cyclical 
period modulation is that it is the result of a solar-type mag- 
netic activity cycle in the secondary star (Applegate & Pat- 
terson 1987; Warner 1988; Bianchini 1990). Richman et al. 
(1994) proposed a model in which the Roche lobe radius of 
the secondary star Rl2 varies in response to changes in the 
distribution of angular momentum inside this star (caused 
by the magnetic activity cycle), leading to a change in the 
orbital separation and, therefore, in the orbital period. As 
a consequence of the change in the Roche lobe radius, the 
mass transfer rate M2 also changes. In this model, the orbital 
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period is the shortest when the secondary star is the most 
oblate (i.e., its outer layers rotate faster), and is the longest 
when the outer layers of the secondary star are rotating the 
slowest. 

The fractional period change AP/P is related to the 
amplitude A(0 — C) and to the period P,-aod of the modu- 
lation by (Applegate 1992), 



(1) 



Using the values of A and C in Table ^, we find AP/P = 
4.4 X 10~^. This fractional period change is comparable to 
those of other short-period CVs, but is one order of mag- 
nitude smaller than those of the CVs above the period gap 
{AP/P ~ 2 X lO-*^) [Warner 1995]. 

The predicted changes in Roche lobe radius and mass 
transfer rate are related to the fractional period change by 
(Richman et al. 1994), 



ARl2 _ 

Rl2 
and by, 

AM2 _ 

Mo " 



39 



1+9 



2/3 



1.22 X 10" 



AO. 

lo-^n 



T'i 



AP 

AO. 



AP 



(2) 



(3) 



where Afl/^1 is the fractional change in the rotation rate of 
the outer shell of the secondary star involved in the cyclical 
exchange of angular momentum, q (= h'h/Mi) is the binary 
mass ratio, and the minus signs were dropped. 

In the framework of the disc instability model (Smak 
1984, Warner 1995 and references therein), the mass trans- 
ferred from the secondary star will accumulate in the outer 
disc regions during the quiescent phase. Hence, changes in 
M2 will mainly affect the luminosity Lt,g of the bright spot 
where the infalling material hits the outer edge of the disc. 
The luminosity of the bright spot is Lbs oc M2. The bright 
spot in Z Cha contributes about 30 per cent of the total 
optical light (on an average over the orbital cycle) [Wood 
et al. 1986]. Therefore, one expects that the changes in mass 
transfer rate lead to a modulation in the quiescent bright- 
ness of the system of Am = (ALts/L) ~ 0.3 (AM2/M2) ^ 
(0.06-0.12) mag for Afl/il = (1-2) x 10"^ (Richman et al. 
1994) and g = 0.15 (Wood et al. 1986). 

Ak, Ozkan & Mattel (2001) reported the detection of 
cyclical modulations of the quiescent magnitudes and out- 
burst intervals of a set of dwarf novae, which they inter- 
preted as the manifestation of a magnetic activity cycle 
in their secondary stars. They found that the quiescent 
brightness of Z Cha is modulated with an amplitude of 
Am = 0.16 mag on a time scale of 14.6 yr. Unfortunately, 
they do not list the epoch of maximum brightness. 

The observed amplitude of the brightness modulation 
is larger than that predicted from the orbital period change 
with the assumption of AQ./il = (1 — 2) x 10~^. The model 
of Richman et al. (1994) can account for both the measured 
period changes and the observed brightness modulation if 
AQ/n ~ 2.7 x 10"^ This yields ARl2/Rl2 1.7 x 10"*. 

The predicted change in the Roche lobe radius of the 
secondary star in Z Cha is comparable to the observed 
change in the radius of the Sun as a consequence of its 
magnetic activity cycle (Gilliland 1981). The period of the 



quiescent brightness modulation is very close to half of the 
period of the observed Pmod and may correspond to the first 
harmonic of a non-sinusoidal or non-strictly periodic period 
change. 

The analysis of Ak et al. (2001) covers only 18 years 
of observations. It would be interesting to see whether the 
28 yr orbital period change also appear as a modulation in 
the quiescent magnitude in a dataset covering a larger time- 
interval. A simple and interesting test of the Richman et al. 
(1994) model is to check its prediction that the maximum 
of the brightness modulation coincides with the minimum 
of the orbital period modulation, which occurred at i5 ~ 
5.2 x 10* cycle (or about JD 2 444 150) for Z Cha. 

Finally, we remark that, if the magnetic activity cy- 
cle explanation is right, our confirmation that Z Cha also 
shows cyclical period changes underscores the conclusion of 
Ak et al. (2001) that even fully convective secondary stars 
possess magnetic activity cycles (and, therefore, magnetic 
fields). 
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